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Abstract 



Variational principles for field theories where variations of fields are restricted 
along a parametrization are considered. In particular, gauge-natural parametrized 
f^ ] variational problems are defined as those in which both the Lagrangian and the 

""p^ ' parametrization are gauge covariant and some further conditions is satisfied in or- 

^ I dcr to formulate a Nother theorem that links horizontal and gauge symmetries to the 

I ' relative conservation laws (generalizing what Fernandez, Garcia and Rodrigo did in 

'^ I some recent papers). The case of vakonomic constraints in field theory is also studied 

C^ ' within the framework of parametrized variational problems, defining and comparing 

two different concepts of criticality of a section, one arising directly from the vako- 
nomic schema, the other making use of an adapted parametrization. The general 
theory is then applied to the case of hydrodynamics of a charged fluid coupled with its 
j^ ' gravitational and electromagnetic field. A variational formulation including conserved 

H I currents and superpotentials is given that turns out to be computationally much easier 

than the standard one. 

1 Introduction 

In some recent literature an interest emerged about variational principles in which some 
specific requirement of physical or mathematical nature force us to restrict the variations 
of fields according to some a priori given rule. 

As an example, let us mention Mechanics with non-integrable constraints on velocities. 
As it is well known there are at least two diff'erent procedures (so called vakonomic and 
non-holonomic) which lead to different equations of motion, it is now clear [18], at least for 
affine constraints, that the more realistic trajectories are solutions of the non-holonomic 
one. However the situation is less clear in more general cases (see [20] and references 
quoted therein), so that we do know very little about constraints in field theory (see 
[23, 17]) and we cannot yet decide whether one procedure is better then the other for a 



reasonable extension out of Mechanics. In any case we shall not enter this question here, 
but just mention that both procedures can be implemented as variational problems with 
constrained variations (see [11]): they differ exactly in the choice of the rule to be imposed 
to admissible variations. 

Another important case which was recently investigated is Euler-Poincare reduction (see 
[2, 3, 9]) where a variational problem with constrained variations arises as a Lagrangian 
reduction of a free variational problem with symmetries. 

A more classical example is presented in the well known book by Hawking and Ellis 
while introducing the equations of motion for relativistic hydrodynamics [13]. There the 
constraint on variations comes from the physical requirement that variations preserve the 
continuity equation. 

Motivated by all these examples, in this paper we first discuss the formal development of 
variational problems in which the variations are parametrized by sections of some vector 
bundle (Section 2). This topic was already addressed in [8] at least for natural constrained 
field theories arising from the vakonomic technique; we recover here, and sometimes gener- 
alize, some of the results of [8] without asking that the parametrization necessarily arises 
from the vakonomic handling of a constraint, and we develop the case of gauge covariant 
field theories (much wider then the natural case) for which we develop a full Nother theory 
to link horizontal and gauge symmetries to the relative conservation laws (Section 3). In 
Section 4 we specialize to the case of vakonomic constraints in field theories and define 
two different concepts of criticality for a section: one arising directly from the vakonomic 
paradigma, the other following from a compatible parametrization. We also show under 
what conditions the two concepts coincide. In Section 5, finally, our results are applied to 
the case of relativistic hydrodynamics of a charged fluid. 

An advantage of the framework we shall propose hereafter is to provide a unifying language 
able to deal with all variational problems with constrained variations, no matter whether 
they come from the vakonomic scheme, from the non-holonomic one, from a reduction 
problem or from any other physical or mathematical requirement. In our hope, therefore, 
this unified language can even allow a better comparison between different schemes. 

2 Variational calculus with parametrized variations 

Notation follows [7], to which we refer the reader for further details. Let C — > M be the 
configuration bundle of a physical theory the sections T(C) of which represent the fields 
(by an abuse of language we will often denote bundles with the same label as their total 
spaces) . 

Definition 1 A Lagrangian of order k on the configuration bundle C — y M [m = dimM) 
is a fibered morphism 

L: J^C — > A'^T*M. 



Definition 2 A parametrization of order s and rank I {s > I \s required) of the set of 



TE 



constrained variations is a couple {E,F), where E \s a vector bundle E — > C while P is a 
fibered morphism (section) 

¥:rC — > (fE)* ®jic VC- 

here J denotes jet prolongation and V is the vertical functor. 

If {x^,y°',e ) are local fibered coordinates on E and {da\ is the induced fiberwise natural 
basis of yC, a parametrization of order 1 and rank 1 associates to a section y°'{x) of C 
and a section e (x,?/) of E the section 

[p\{x''J{x),d,y\x))e''{x,y{x)) + p\ ^{x^ J{x),d,y\x))d^e^{x,y{x))\ da 

oiVC. 

Definition 3 Given a compact submanifold D C M with suitably regular boundary dD, 
an admissible variation (for a Lagrangian with order k and a parametrization of order s 
and rank /) of a configuration a G r(C) on I? is a smooth 1-parameter family of sections 
{^t}te\^i,i[ C T{tt-^D) such that 

1. do = ctId 

2. vtG]-l,l[, /-Vi|afl=/-V|aD 



3. there exists a section e G T{E) such that J^(Tt|t=o =< P|i e > oj^cr and j'^+' ^(e 
^)|OD = 0. 



o 



Definition 4 The set {C, L,P} where C — > M is a configuration bundle, L is a Lagrangian 
on C and P is a parametrization of the set of constrained variations is called a "parametrized 
variational problem". 

Definition 5 We define critical for the parametrized variational problem {C, L,P} those 
sections of C such that, for any compact domain D C M and for any admissible variation 
{at} defined on D one has 



= 0. 

i=0 



Accordingly, if we use the trivial parametrization P : C — > VC* (8)c ^C that to any p £ C 
associates the identity of VpC then the third condition becomes empty and we recover free 
variational calculus. 

For an ordinary variational problem with Lagrangian L = Cds, criticality of a section of C 
is equivalent in local fibered coordinates {x^,y'^) to the fact that for any compact D d M 
and for any V G T{VC) such that j^~^{V o (t)!^^) = one has 



D 



dC dC dC 

Q^V\x, y{x)) + Q-d^V\x, y{x)) + ■■■ + ^- d^, ■ ■ ■ df,,V'{x, y{x)) 



ds = 0. 



Explicit calculations (see [7]) show that the local coordinate expressions given above glue 
together giving rise to the global expression 

yO C Af, Vy G VC s. t. 3^-^{V o a)\QD =0, [ <SL\ j^V > oj^a = (1) 

JD 

where 5L is a global fibered (variational) morphism 

5L : J^C -^ {J^VCy ®jkc A"'T*M. (2) 

The rigorous definition and a short account on variational morphisms will be given in 
Appendix A, where the reader can find some non-conventional definitions that we have 
generalized in order to embrace the case of parametrized variations. 

2.1 P-first variation formula and P-Euler-Lagrange equations 

To define criticality for parametrized variational problems with a /c -order Lagrangian 
and a paranietrization of rank /, we have to restrict variations to those V G VC in (1) 
with J (V^ o (7)|ai5 = that can be obtained through the paranietrization from a section 
e oi E satisfying ^"'"'^"^(e o a)\gj;) = 0. Let us make t explicit on a first order Lagrangian 
and a paranietrization with rank 1 and order 1: if (x^, y"", e ) are local fibered coordinates 
on E and {p'X^^ ~^ Pa ^'^A«^'^)'9a is the local representation of < P \j^e > oj^a, criticality 
holds if and only if for any compact D C M, for any section e with coordinate expression 
e {x,y) such that both e {x,y{x)) = and d^e (x,y(a;)) = for all x G dD, we have 



L [^^^^'^ ^^^ '"^"'^^ + ^MPaS^+Pa ■'d.e^) 



ds = 0. (3) 



To set up an intrinsic characterization of critical sections possibly given by a set of differ- 
ential equations let us introduce the following procedure. 
Let us take any paranietrization P and think of it as a morphism 

P' : rC X fE — > VC 

linear in its second argument. Take its /c -order (holonomic) jet prolongation (see Ap- 
pendix A, Definition 35) 

j^r : r+^c X f+^E — > J^VC 

and read it as a bundle morphism 

j¥' : r+^C -^ {f+^Ef ®js+kc J^VC 

(for simplicity we still denote it by j P', but the same abuse of notation cannot safely be 
adopted also for P and P' because the jet prolongations of the two objects give completely 
different results). 

Let us consider the variational morphism 5L : J C — > {J V C)* (^ jk(j K"^T* M introduced 
in formula (2) and take its formal contraction (see Definition 37) with j^¥' . What we get 
is the variational i?-morphism 



Criticality of a section a can now be recasted in the following global requirement 

VD C M, Ve G r(^) s. t. j^+''~^{e o a)|g^ = 0, 



< 6L I j'^F' > 

D 



jl + k\ ^ jk + S^ ^ Q 



This intrinsic procedure leads in the case of first order Lagrangians and first order first 

rank parameterizations to the local expression (3). 

Now we can try to turn this requirement into a set of differential equations. The first step 

is to search for a first variation formula for parametrized variational problems. 

Thanks to Theorem 45, for any choice of a fibered connection (7, F) on E — > C — > M 

(see Definition 39 and 41), there exists a unique pair (IE, F) of variational morphisms 



E(L,P) : j2fc+'+^c — > E* (^c ^"'T*M 

F(L,P,7) : J^k+l+s-l^ _^ (/+^-l^)* ®Jl+,-^c^ 

reduced with respect to (7, F) (see Definition 44) such that Ve G F(E') one has 

<5L|j¥'> y+'=e\ =<E|e> + Div<F|/+''-^e> . (4) 

In the previous equation Div stands for the divergence of a variational morphism (Ap- 
pendix A, Definition 38) that coincides when computed on a section with the exterior 
differential of forms. 

Definition 6 We call F-Euler- Lagrange morphism the variational morphism E,{L,S,F) and 
F-Poincare -Cartan morphism the variational morphism F(L,5', P, 7). 

Proposition 7 A section cr of C is critical for the parametrized variational problem (C, L,P) 
(with L of order k and P of order s and rank /) if and only if < E | e > of^+^+^a = 0. 



■l + kA ^ -k + S^ ^ Q 



Proof: We have 
VD C Af, Ve G T{E) s. t. i^^^~^{e o (t)\q^ = 0, / / < 5L\j^F' > 

t 
yD C M, Ve G T{E) s. t. f+''-\e o a)|g^ = 0, 

/ < E I e > oj2'=+'+V + [ d{<¥\ j'+'^-ie > of^+'+^~^a) = 
Jd Jd 

t 

E o j2fc+«+s^ ^ 

The first equivalence is ensured by the first variation formula (4) and the Definition 38, 
while the last holds in force of the fact that j^~^^~^(e o o") L„ = and that e and D are 
arbitrary. Jl» 



We stress that the number of derivatives of fields in E and F is not sharp: in fact, due 
to the particular structure of the morphism < 6L\j^¥' >, it is actually always lesser 
then the one expected from the naive issue. Moreover, the morphism E(L,P) in the 
splitting is truly unique (i.e. independent on the connection used), while F(L,P, 7) is just 
one representative in the class of variational morphisms fitting into formula (4) chosen o 
that it is reduced with respect to the fibered connection (7, F); more precisely, the final 
expression of the morphism F(L, P, 7) depends at most on the spacetime connection 7, 
while ifl + k < 2 (the most common case) also F can be defined independently on any part 
of the fibered connection, so that the splitting algorithm can be performed in coordinates 
simply by integration by parts on the derivatives of the sections of E. 
In the case k = 1, 1 = 1 and s = 1, with the same coordinate representation as in formula 
(3), the morphisms E and F are the following: 



< E I e > of a 



< F I fe > of a 






Qya ''''dy-l^^ 



e^ ds 
ds 



and they can easily be obtained integrating by parts the integrand of equation (3). 

2.2 Symmetries 

The second natural question arising in our investigations is whether the link between 
symmetries and conserved currents given in the standard Nother theory is preserved if we 
vary configurations only along the parametrized variations. First of all the definition of 
symmetry for a parametrized variational problem has to be changed with respect to the 
ordinary one in order to ensure the invariance of the parametrization as well as that of 
the Lagrangian. 

Definition 8 Let C ^^ M be a configuration bundle, L a Lagrangian of order k and {E,¥) 

a parametrization of order r and rank / on the vector bundle E — ^ C. 

Let Z : E — > TE be a vector field on the bundle E — > M such that Y = T-keZ is a vector 

field on C , while X = TtteY is a vector field on M. 

Calling ^f 3, "^Ys ^'^^ ^y? ^^^ flows of fibered transformations generated by the vector fields 

Z and Y on the bundles E, C and VC respectively, we can drag any section e € T{E) 

along ^zs getting the section Sg = (^f g)*e, any section a G F(C) along ^y^ getting 

as = {^Ys)*^ ^''^^' ^°'' ^i^y section V G T{VC), we can also drag V o a along ^y^ getting 

i^^^rvoa. 

The vector field Z is an infinitesimal Lagrangian symmetry of the (parametrized) variational 
problem {C, L,P} if the following covariance identity 

<5L\f£Y(J>= D\M{ixL) + D\y < a\fY > (5) 

holds for some morphism a : J^+'^C — > {J^TC)* ®jkc K^-^T*M, and moreover we have 

< P I j'e, > of as = (^1^,?)^(< P I j'e > of a). 



The definition of infinitesimal symmetry lias now been extended to variational problems 
with parametrized variations. Nevertheless in order to establish a correspondence between 
symmetries and conserved currents we need some more restrictive hypothesis. Usually (in 
the case with free variations) one considers the covariance identity (5) and simply splits 
6L in its Euler-Lagrange part plus the divergence of the Poincare-Cartan part: being the 
first vanishing on shell, one gets the conserved current as the argument of the residual 
divergence. In our framework, on the contrary, field equations arise from the splitting 
oi < 6L\ j^F' > so that we would need a rule to associate to any symmetry Y and to 
any configuration a a section of the bundle of parameters that in turn is mapped by 
the parametrization into £yO'. Without any further hypothesis such a section is not 
guaranteed to exist. In thee next Section we will then state all the necessary assumptions 
which are needed to get a Nother theorem in the case of gauge-natural symmetries of 
gauge-natural variational problems in the next section. 

3 Gauge-natural constrained theories 

Gauge-natural field theories have been shown to be the most general setting to describe 
Lagrangian field theories with gauge symmetries, including gauge theories. General Rela- 
tivity in its different formulations, bosons, spinors and also supersymmetries (see [5, 7]). 
In the book [7] it was also described in detail how to cope with conservation laws linked 
to pure gauge (vertical) symmetries and how to implement the conservation of energy 
and momentum as Nother currents relative to an horizontal lift of space-time diffeomor- 
phisms. The goal of this Section is to extend the previously known results to the case 
of constrained variational problems under suitable conditions. The original material on 
gauge-natural bundles can be found in [4]; a standard reference with many theoretical 
improvements is [16]; a more applicative approach can be found in [7], while a friendly 
introduction to the subject is also provided in Appendix B. 

3.1 Morphisms between gauge-natural bundles 

In this sub-Section we recall the fundamental definitions about morphisms between gauge- 
natural bundles, their push-forward along gauge transformations and their Lie derivatives. 
Gauge-natural morphisms are also defined. 

For any principal automorphism ^ = ((/>, ip) of the principal bundle P — *• M let us denote 
by !B(^) = {4>, ^^b) the gauge-natural lift of ^ to the gauge natural bundle ^{P). 

Definition 9 A bundle morphism between gauge-natural bundles M : 53(P) — > 2)(P) 
projecting onto the identity is said to be gauge-natural if for any local principal automorphism 
^ of P the following diagram commutes 



^{P)^^^{P) 



<B(*) 



S)(^) 



^{P)^^T>{P) 



Definition 10 Let M : 53(P) — > ^{P) be a vertical (i.e. projecting onto the identity) 
morphism between gauge-natural bundles. We define \ts push-forward along the local principal 
automorphism ^ = ((/>, V') to be the unique vertical morphism ^*M : 53(P) — > 2)(P) such 
that for all section a G r(*B(P)) we have 

2)(*)*(M o cr) = ^*M o «B(^)*o-. 

As a consequence we have the fonowing exphcit rule of calculation: 

^"^M = D(^) o M o 53(1'-^) = vf J, o M o ^-^ 

Definition 11 Let M : 53(P) — > 2?(^) be a vertical morphism between gauge-natural 
bundles over the same base and ^^ = {(ps,ijjs) a 1-parameter family of local principal auto- 
morphisms having E € r(IGA(P)) as generator. We define the Lie derivative of the morphism 
M along E to be the morphism £eM : «B(P) — > VD{P) that fulfills 



ds 



s=0 



From the definition it follows immediately that the expression of the Lie derivative of a 
morphism is linked to that of the Lie derivative of sections by means of the following rule: 
for ah a G r(«B(P)) 

I^M oa = £3{Moa)-TMo £^a. 

If we use coordinates (x^, y") on 53(P) and {x^, z ) on S{P), the local expression of M is 

Let us call Ssg and Hj) the gauge-natural lifts of an infinitesimal generator of principal 
automorphisms E respectively to ^{P) and Ti{P), and let 

A _ „^ <9 Aa "9 A _ „u 5 AA f^ 

be their coordinate expressions. 

The Lie derivative of M has thence local coordinate expression 



■^2) - ir:^-'B 



Remark 12 Let ^(-P) and Ti{P) be gauge-natural bundles of order (r, s) and (j, A;) respec- 
tively. We can think at the Lie derivative as a fibered morphism 

£m : 'B(P) -^ (J'"IGA(P))* (g)j,(P) VDiP) 

with m = maxj?', j} such that if H is vertical (i.e. a section of ^ ^-> IGA(P)), £-=yi depends 
in fact on the derivatives of E only up to the order / = max{s, k}. 



Proposition 13 A necessary and sufficient condition for a vertical morphism between gauge- 
natural bundles M : 5S(P) — > ^{P) to be gauge-natural is that for any 1-parameter family 
of local automorphisms ^^ of P generated by H, one has 

^*M = M or £sM = 0. 

Proposition 14 Any global morphism M : S(P) — > 2)(P) has to be gauge-natural and 
any local gauge-natural morphism can be extended to a global one. 

Proof: (Sketch) A trivialization of P induces both a trivialization of 53 (P) and 2)(P) 
(see [7]) and gauge-naturahty is exactly the same as invariance with respect to changes of 
triviahzations of P. X 

3.2 Variationally gauge-natural morphisms 

In this sub-Section we present some technical properties that follow from requiring a local 
Lagrangian not exactly to be a gauge-natural morphism (that would imply globality) 
but to satisfy a slightly more relaxed condition that amounts to say that every gauge 
transformation has to be a Lagrangian symmetry according to definition 8 if a gauge- 
natural parametrization is provided. 

This generalization is needed in order to deal with the case presented as an example in 
Section 5.1. 

The Theorems we shall develop here can be applied to ordinary variational calculus, but 
they have been introduced with the explicit aim of embracing the case of parametrized 
variations; their consequences will be analyzed in the next sub-Section. 

Definition 15 Let £(-P) — > M be a gauge-natural bundle and F — > ^{P) be a vector 
bundle such that the composition F — > ^{P) — > M is a gauge natural bundle on M that 
we call ^{P). A local variational morphism 

M : j''it{P) — > {J^d{P)T ®j<^c A™""r*(M) 

is said to be variationally gauge-natural if for any 1-parameter family of local automorphisms 
^s of P generated by H there exists a 1-parameter family of variational morphisms {««} with 

Us : j'^-^diP) — > (j'^-i5'(-P))*®A™""~^r*(M) 

such that, for every section X G T{^{P)) one has 

** < M I j^X >=< M t ^^i'^X > -h Div < a^ I ^*j^~^X > . 



Remark 16 In terms of Lie derivatives there are many possible equivalent infinitesimal for- 
mulations of the definition of variational gauge-naturality. Let us consider the morphism 

M : J^(t{P) — > iJ^d{P)T 0jhc A'"""r*(M) 



and let us think of it as a morphism 

M' : j''<t{P) X j'^diP) — > A™-"T*(M) 

linear in its second argument. Let us suppose that C(P) is gauge-natural of order {r,s) and 
er(P) of order {j,k) and that m = min{r, j}, according to Remark 12; then the Lie derivative 
of M' with respect to a gauge generator H is a morphism 

£eM' : j'=e:(P) X j'^diP) — > (J™IGA(P))* 0m FA™~"r*(M), 

linear in the second argument, that can also be thought as 

£eM' : J^(t{P) — > (j'^diP))* ® (J™IGA(P))* ® yA™-"r*(Af) 

where the tensor product is on the base J^(t{P) (for simplicity we still denote it by £^M'). 
Using the map tt2 introduced in Remark 54, we can say that variational gauge-naturality of M 
is equivalent to the existence of a morphism a : J^~^<t{P) — > {J^^{P))* ® (J™IGA(P))* 
VK^-^T*{M) such that 



f~'x 



TT2 < £eM' I fX >= Div (vrs/ < a I j'^-^E > 
and the link between a and the 1-parameter family {as} of definition 15 is 

d 



ds 



as\s=Q=<a \j'" ^=> . 



Concrete examples in local coordinates of how to check variational gauge-naturality of a 
Lagrangian morphism as well as that of an Euler-Lagrange morphism will be given in the 
proof of Proposition 19. 

In our exposition, two technical lemmas precede the main basic properties of variationally 
gauge-natural morphisms and in particular of variationally gauge-natural Lagrangians in 
order to simplify their proofs. The first lemma is very well known, but we report it here 
just to display some coordinate expression which will be easily recognized later. To our 
knowledge the other statements have not yet been proved in the same generality. 

Lemma 17 If L = Div A then for every X G T{VC) there exists a morphism 5K such that 
< 5L I /X >= Div < M I f~^X >. If A is global, so is 5A. 

Proof: Let A be a local variational morphism 

where M is the base of the configuration bundle C. Let us use coordinates {x^, y%) on J^C 
(;U runs from 1 to dim M, while a is a multiindex of length < |a| < k) and let us express 
A in coordinates as A = A'^(x'^, y%) ds^ (with < |q;| < k) while L = d^A.^{x^^ ,y%) ds. 
The coordinate expression for the morphism < 6L\j^X > is the following 

8d ^^^ 
<5L\j'X>=^^XBds. (6) 

dy$ 

10 



where the multiindex v has length < |i/| < fc. 
We have 

where the multi-index a has length \ol\ = \v\ — 1 and by 1^ we mean a multiindex with 
1 entry in position // and zero elsewhere, while for |i>| = there is no q to fit into the 

whil 



formula so that one has to set by ** = 0, while for \v\ = k we simply have 



'^M = '■ 



Substituting now into (6) we get 



<bL\fX>-- 






We have found the local expression of the morphism 5A 



n ds = d,\ -^Xl 1 ds. 



bk = {bkr ds^ with {bkY = -^^X% 



dkt" 

of which we wanted to prove the existence. It remains to straightforwardly check that if 
A is global so is bk; and that in this case the previous one is an invariant expression that 
does not depend on the chosen coordinates. X 

Lemma 18 The gauge-natural lift j'^'^va: of 3" infinitesimal generator H G r(IGA(P)) of 
principal automorphisms of P to the A:*'^-jet prolongation of the vertical bundle V^{P) of a 
gauge-natural bundle C(P) can be computed from the lift j^Hg; to the same jet prolongation 
of the base bundle £(-P) according to the following rule: let (a;^,yp) be fibered coordinates 
on J''€{P) and (x^,^^,^) on j'=Fe:(P) [u is a multi-index of length < |i>| < k). If 






is the coordinate expression of j Hg then the coordinate expression of j Eyir is 

' -""^-^ d^^^^' -^^'M ^' ^^^^'M 

with 

where j2 \s a multi-index of length < |/u| < k. 

Proof: If we have a gauge-natural bundle ®(-P) with a system of local coordinates 
(x^,p ), given the gauge-natural lift H(g of any infinitesimal generator H of principal 
automorphisms of P with local coordinate representation 






11 



then the hft Hytg on the vertical bundle of *B, where we use coordinates {x'^,p^,e^), has 
the following local representation 



^x^^ dp^ dp^ de^ 

Now plug in <B(P) = J''C{P) and, thanks to the isomorphism ViJ'^aiP)) ^ j''{V<i{P)), 
one gets the thesis. Jl» 

Proposition 19 If L is a variationally gauge-natural (local) variational morphism then the 
gauge-natural morphism 6L introduced in formula (2) is variationally gauge-natural, too. 

Proof: Let us introduce a fibered coordinate system {x^,y"') on the gauge-natural 
configuration bundle ^(-P). On J (t{P) we can naturally induce the system of fibered 
coordinates {x^,yp) where i/ is a multi- index of length < |z?| < /c. 
For a A;-th order Lagrangian and for any section X of the vertical bundle VC{P) we have 



dy$ 

Let us adopt the same notation of the previous Lemma for the lifts of infinitesimal gener- 
ators of automorphisms of P to the configuration bundle and its vertical bundle. 
Furthermore let (r, s) be the order of C, so that the order of J^'^L is at most {r + k,s + k). 
The Lagrangian L : J €. — > N^T*M is variationally gauge-natural so that a morphism 
a : J^-^(t{P) — > {r+^-HGk{P)y VA'"-"r*(M) exists such that 

TTzi^HM' = Div (1:2 <a\ j'^-^H > ) = Div A. 



having called A the morphism A = 712 < a | j'^+^'^^H >; let A'^ be its components with 
respect to the decomposition A = A'^ds^. 
We have in coordinates: 

which is equivalent to 

Thanks to Lemma 18 we can notice that if the order of (t(-P) is (r, s) then the order of 
VdiiP) is {r + l,s + l). 

Gauge-naturality of 6L : J^^{P) — > {J^V€{P))* A"'T*M means that there exists a 
morphism /3 : J^-^(t{P) — > {J^^Wt)* ® {J'-+HGA{P))* ® FA"^-"T*(M) such that 

vra < i^sM' | j^X >= Div (vra^ < /3 | f+''E > j^~^X 

Let us then compute < i^^M' | j^X >. We have (both multiindices v and ]1 have length 

< |z^| = |/i| < k). 



12 



r) BT 

vra < £^M I j^X >= <j X?,-^—-d^£sy' - d. 



-n 



d dL 



d dL 



(S^''- 



+ 



(j'Sc) 



dL 






BT 

+ ^-^da^X^ \ ds 

dL 



dy'^pdx^ dy^^dyl 

Integrating by parts the first two addenda and applying (7) we get 



d^[Xl~''\+^_da=.X^\ ds. 



^, < £sM' I j'X >= \ X? A [Lrf^H- + d„A-] + ^^^_U'^e)l+ 



-^x^d^- 

dyl ' 



—^S'd^xl + :^d^£^x' \ ds 



V + -;r^dy£-=X 
dyl 

dL ,,„ d , .u- sb dL , .i, - ,6 ,,„ d 



dy. 



n^aij'^^). - 7i3(/^v^c); + X?— (d.A-) \ ds. 



dyl ^dy 



dyl 



'dyl 



Thanks to Lemma 18 the first two addenda cancel, while thanks to Lemma 17 the argument 
of the derivative in the third addendum is the divergence of the morphism 5 A. with local 
expression 



< 5K I j'~'X >= ( Xf^ ) ds. 



dy't. 



Thus we have 



vr2( < /3 I Z+^^H > j'^-^Xj =< 6A \ j^-^X > 
and 5L is variationally gauge-natural. 



Theorem 20 The volume part of a variationally gauge-natural (local) variational morphism 
is gauge-natural (and thence global). 

Proof: Let P — > M be a principal bundle, G^(-P) a gauge-natural bundle and ^iP) a 
gauge-natural vector bundle. Let M be a gauge-natural morphism 

M : j''(t{P) — > {J^d{P)y 0jhc A™T*M 

that for ah V € T{^{P)) splits into 



< M [ j^V >=< V I y > + Div < B I j'^-'^V > . 



with 



V = V(M) : f'+^^{P) — > [^{P)T ®^K^T*M 

B = B(M,7) : J^+^-^<t{P) — > {.J^-^^{P)Y (S)jh-i^N^-^T*M. 

For any 1-parameter family of local automorphism ^ s of P and for all V S T{'^{P)) we 
have 



'f * < M I /y >= ** < V I y > +^^ 



Div < B 1 j'^-^V > 



As it happens whenever we consider a natural bundle on M as gauge-natural on any 
P — > M the gauge-natural lift on K^T*M of any gauge transformation coincides with 
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the natural lift of the diffeomorphism onto which the gauge transformation projects, so that 
if ^s = {(ps,ips) then A'"r*(^s) = A'^T*{(ps)- If computed on a section, the divergence 
coincides with the exterior differential of forms (see Definition 38). Accordingly, as the 
push-forward along a diffeomorphism commutes with the exterior differential we have the 
following splitting 



^* < M I fV >=< ^*V I ^^y > + Div 



^^ < B I j^- V > 



By variational gauge-naturality of M we also have a 1-parameter family of morphisms 
{os} such that 

^* < M I fV > =< M I ^^/y > + Div < a, I ^*(i^~V) > 

=< V I ^*y > + Div < B + a, I "^tif'^y) > ■ 

Thanks to the uniqueness of the volume part of the splitting of a variational morphism 
we can conclude 

that ensures gauge-naturality of V. X 

3.3 Gauge-natural variational calculus with parametrized variations 

Definition 21 A gauge-natural variational problem with parametrized variations is defined 
by the set {<t{P), L,^,FJ,u;) of the following objects: 

1. a gauge-natural configuration bundle ^{P) — > M of order {r,s) with structure G- 
bundle P -^ M; 

2. a variationally gauge-natural (local) Lagrangian morphism L of order k; 

3. a vector bundle F — ^ ^{P) such that the composite projection F — ^ M is a gauge- 
natural bundle d{P) of order {p, f) called bundle of parameters; 

4. a gauge-natural parameterizing morphism P : j'^(t{P) — > {Ji^{P))* ®jqe{p) V(S:{P) 

5. a gauge-natural morphism J : J^d — > (J'^~«IGA(P))* (E)jic5(-P) such that we have 

(P I f < Jlf-^H > > o jV = £^a; 

6. a morphism w : J^ft{P) — > C{M) x^/ C{P) that associates to any configuration a 
couple (7,0;) where 7 is a linear connection on M and w is a principal connection on P. 

Let us remark that the connection 7 is needed to construct the P-Poincare -Cartan mor- 
phism F(L,S', P, 7) only \i k + I > 2, otherwise this requirement can be dropped, while 
the connection uj on P will be needed in every case to distinguish horizontal and vertical 
gauge symmetries (as we will see in a while). 
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Theorem 22 A gauge-natural variational problem with parametrized variations {<t{P),L, 
5, P, J) leads to gauge-natural (global) Euler-Lagrange equations. 



Proof: According to the splitting 



j^+^e) =< E I e > + Div < F I j'+''~^e >, 



< 5L I j^¥ > 

the Euler-Lagrange equations of a gauge-natural variational problem with parametrized 
variations arise as the boundary part of a variationally gauge-natural variational mor- 
phism; in fact 5L is variationally gauge-natural by Property 19, and so is P by hypothesis 
(hence also its jet prolongation). Their contraction is thence gauge-natural by definition. 
Thanks to Property 20 we can conclude that the equations are gauge-natural. X 



Let us now remark that according to the fact that the Lagrangian is variationally gauge- 
natural, every gauge-natural lift on ^{P) of an infinitesimal generator H G r(IGA(P)) of 
automorphisms of P (or, in a more physical language, every gauge transformation) turns 
out to be a symmetry of our variational problem. With the help of the other requirements 
of Definition 21 we are able to associate to any of these gauge symmetries a gauge invariant 
"Nother" conserved current according to the the rule that follows. 

Theorem 23 -Nother theorem- Let us consider a gauge natural variational problem with 
parametrized variations ((J^(P), L,5, P, J) (orders are as in Definition 21) such that 

■K2£=.L = Div (vrs < a I f+^-^r. > 

The unique fibered morphism (momentum map) 

£ : j2'=+'+^-ic:(P) — > J^-«^i|GA(P) ® A"-^T*M 
such that for all a € T{it{P)) and for all E G r(IGA(P)) one has 

< £ I f~i~^E > of^+^+'-^a = 

< F I y+^-ij oE>-TpE^L-7r2<a\ f+'^-^E > ) o f^+i+s-i^ 



associates to any couple {a, E) formed by a section a G T{(t{P)) and an infinitesimal generator 
of principal automorphisms E € r(IGA(P)) a current (i. e., a (n — l)-form) that is conserved 
(closed) on-shell. Moreover the morphism £ is gauge-natural (global). 

Proof: By gauge-naturality of L : J C — > A™T*M one can find a local morphism 
A : J*'- 1(2: — > yA'"-ir*M such that 

tt2£eL = Div (vra < a I f+''-^E >). 
7r2£EL =< 6L\j''£e(t > - Div (TpHjL) = 



■l+k TT I -r-q 



< SL j f¥' > f^" < J I f-'^E > ) - Div {TpEjL) 



E 



< J I f-'^E > ) + Div( F 
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i'+^'~i < J I f-'^E >) - Div (TpEjL) 



Thence 



< J 1 /"^ 


^>) 


= 


= DIv/f 


j'+^-i < J 1 f-'^E > 




-(E 


< J 1 f'E > 



- Div (TpEjL) - Div (vra < a | /'+'=- ^H >) 

so 

= Bw<£\ f-'^-^E > . (8) 



While from the definition of £ one can hardly argue how many derivatives of E will appear 
in its final expression (cancellations may arise!) from the latter identity (that holds off- 
shell) we can say that E can enter with its derivatives up to the order r — q — 1. If o" is a 
solution of the corresponding Euler-Lagrange equations then we have 



d 



F 



'+'^~^ < J j f-'^E >)- TpEjL - vra < a I f+''~^E > ) o j2fc+«+s-i 



0. 



To prove gauge-naturality of £ let us notice that identity (8) holds; thence, as the push- 
forward and the exterior differential of forms do commute, it is sufficient to prove gauge- 
naturality of the left hand side. The morphism E is gauge-natural by Property 20, and it 
is contracted on J that is gauge-natural by hypothesis; thus the thesis follows. X 

In the natural case we were able to construct Nother currents related to diffeomorphisms, 
and their integrals on (m — l)-surfaces were linked to the energy-momentum of the system. 
In the gauge-natural case we cannot any longer lift diffeomorphisms of M, but according to 
[7], if on P we have a connection gauge-naturally derived from fields, for any infinitesimal 
generator of principal automorphisms we can define the horizontal part with respect to 
this connection, and compute the conserved Nother currents relative to it. This should be 
regarded as the density of energy momentum in a gauge-natural variational principle also 
in presence of a parametrization relative to a constraint. 

4 Constrained variational calculus and vakonomic field the- 
ories 

Let us now apply the previously developed formalism to field theories with constraints 
within the vakonomic approach. The non-holonomic technique can be studied within the 
same formalism too, but with a different characterization of admissible variations; this 
case will be considered in a forthcoming paper. 

4.1 Vak-criticality 

The definitions we will use to implement the principle of constrained least action are a 
generalization of the ones given by Arnold and others in [1] in the case of Vakonomic 
Mechanics. We will abandon their line when they introduce Lagrange multipliers to char- 
acterize critical sections. 
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Definition 24 Let C -—> M be the configuration bundle. A constraint of order p and 
codimension r is a submanifold S C J^C of codimension r that by tt^^i projects onto the 
whole JP-^C. 

Definition 25 A configuration a € r(C) is said to be admissible with respect to S if its 
p'^'-order jet prolongation lies in S. The space of admissible configuration with respect to S is 

TsiC) = {aeTiC)/\mifa)eS}. 

cpt 

Definition 26 Given a compact submanifold D C M, a vak- admissible variation (at order 
k) of an admissible configuration a G r5'(C) is a smooth 1-parameter family of sections 
We}ee]-iA[ C T{tt''^D) such that 

1. do = a\D 

2. VeG]-l,l[, /-V,|a,5=/-V|a^ 

3. ifae\e=oeJWCnTS. 

Definition 27 The set {C, L, 5} where C — > M is a configuration bundle, L is a Lagrangian 
on C and 5 is a constraint is called a "constrained variational problem". 

Definition 28 We say that an admissible section a € ^s{C) is vakonomically critical (or 

cpt 

vak-critical) for the variational problem {C,L,S} if VI? C M for any vak-admissible variation 
{ae}ee]-i,i[ e T{7r-^D) we have 



d > ^ .fc. 



de Jd 



= 0. 

e=0 



An equivalent infinitesimal condition is that 

VD C M, yV G VC s. t. /y G TS and i'^" V|az) =0, f <SL\ j^V > of a = 

Jd 

where SL is the fibered (variational) morphism defined in (2). 

4.2 Parametrized variations and vakonomic criticality 

In order to turn the vak-criticality condition into a differential equation, the classical 
strategy (at least in Mechanics) is to introduce Lagrange multipliers (see [1]). Here we 
try to profitably use a suitable parametrization of the set of constrained variations. This 
idea was first formally developed by Fernandez, Garcia and Rodrigo in the very interesting 
paper [8]. 
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Definition 29 Let C -—> M be the configuration bundle. A parametrization 

¥s : J'C — > (fE)* ^jsc VC 

of order s and rank / (s > Z is required) of the set of constrained variations is said to be 

vakonomically adapted to the constraint S C J^C if for all e G ^{E) and a G ^s{C) the 
vertical vector field f{< P | j'e > o/o-) e TS*. 

Definition 30 A parametrization P5 : J^C — > (j'E')*(8)jacl/C vak-adapted to a constraint 
S is said to be faithful on a e Ts{C) to S if for all V eVC such that both jP{V oa) eTS 
and j''~^iV oa)\gD = hold, there exist a section e G r(£^) such that < ¥\fe > of a = V 
zn6f-^+\eoa)\9D = ^ 

The fundamental problem of the existence of a (possibly faithful) parametrization vako- 
nomically adapted to a constraint S has not yet been studied in general and will not be 
faced here. For the moment we limit ourselves to consider a given specific parametrization 
as a part of the variational problem, deferring the study of the general case to further 
investigations. The same was done in all the previous papers on this topic (see e.g. [8]). 

Once we have a vak-adapted parametrization P5 we can study P^-criticality, constructing, 
as we have shown in Section 2.1, the relevant P^-Euler-Lagrange equations. 
It is important to remark that even having defined the Lagrangian on the whole configura- 
tion bundle, the first variation formula arising from it does not depend on the value of the 
Lagrangian outside the constraint. In fact if one has a Lagrangian L : J C — > A*" T*M 
such that L = L + N with N proportional to the equation of the constraint 5, then on the 
constraint 5 it automatically holds < 5N \ j^V > oj^a = for all V G VC s.t. fV G TS 
and for all a G Ts{C) therefore the addendum N does not contribute to the first variation 
formula, neither in the field equation nor in the boundary part. 
The link between vak-criticality and P5-criticality is given by the following Proposition. 

Proposition 31 Let a G ^s{C) be a vakonomically critical section for the constrained 
variational problem {C,L,S}; then for any adapted parametrization P^ a is Pg-critical. 
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Proof: We have 



a G Ts{C) is critical 
t 



1 f' 

yo'^C M, V adm. var. {as}, — / Loj^cFe 

de Jd 



= 

£=0 



t 



VD C M, yV e VC s. t. fV G TS and j^~^(j\dD = 0, [ <SL\ j^V > o/ct = 

Jd 



l+kA „ -k+s 



\/D C M, Ve G r(S) s. t. 3^^^~^e\QD = 0, f /<6L\ /P' > 



j'+'^e ) o j^^+V = 



t 

Step (a) is not an equivalence because there can be admissible infinitesimal variations 
vanishing at the boundary with their derivatives up to the desired order that do not come 
from sections of the bundle of parameters that do vanish on the boundary. The last 
equivalence holds in force of Stokes's theorem, the vanishing of j ~^ e on the boundary 
and the independence of the generators of E. J(t 

Corollary 32 Let a G ^s{C) be an admissible P^-critical section for the constrained varia- 
tional problem {C, L, S} and let also Pg be faithful to S on a; then a is vakonomically critical 
for the constrained variational problem {C,L,S}. 

5 Examples: 

To our knowledge, parametrized variational problems were first formally investigated by 
Fer-nandez, Garcia and Rodrigo [8] in relation with their application to Lagrangian reduc- 
tion and in particular to Euler-Poincare reduction [2, 9, 3]. They were already implicitly 
used in the literature to deal with Vakonomic Mechanics [1, 11] and Relativistic Hydrody- 
namics [13]. What we want to present here is a gauge-natural example that also helps to 
clarify the distinction between a vak-critical section of a constrained variational problem 
and a P^-critical section, of which we study the equations of motions and the conserved 
Nother currents. 

5.1 A gauge-natural example: charged general relativistic fluids. 

A charged general relativistic fluid can be thought as a congruence of curves filling a causal 
domain D C M, each of them carrying some rest mass and charge, being source of gravity 
and of electromagnetic interaction. 

Our kinematical description is inspired by [13], where the gravitational field is represented 
by a Lorentzian metric with components g^,y, the electromagnetic potential by a connec- 
tion Afj, on a principal C/(l)-bundle P and the fluid degrees of freedom are represented by 
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a nowhere vanishing 3-forin J = J^ds^ such that the unit tiniehke vector u'^ = J'^/j J| is 
tangent to the flow hues and that the volume form y\g\p = \J\ represents the matter den- 
sity in D. For what concerns the charge density we simply assume that at each spacetime 
point it is proportional to the matter density by an "elementary charge" q. 
An alternative (more common) approach to describe the fluid degrees of freedom (see 
[22, 14, 21]) is to associate to any space-time point x^ three scalar flelds R'^{x^) (a = 1, 2, 3) 
identifying the "abstract fluid particles" . We remark that using the tangent vector to flow 
lines as fluid variable, we loose the direct kinematical identification between different points 
of the same flow line. In fact, having a conflguration J^{x), in order to know whether two 
points of spacetime are connected by the same flow line, we have to solve the ODE that 
defines the integral curves of u^ with initial condition in one of them. 
In order to implement conservation of matter we only allow for closed J, so that for any 
closed 3-surface S in D the fiow of J through E is zero. The bundle of configurations is 
thence 

e:(P) = Lor{M) xm C{P) xm A^T*M 

that is gauge-natural having P — > M with fiber U{1) as structure bundle. The field J, 
moreover, is subject to the constraint S C j^A^T*M represented by the equation dJ = 0, 
or locally d^J^ = (we also require J to be time- like with respect to the unknown metric, 
but this just selects an open set in the configuration bundle, so that the condition is 
automatically preserved by any infinitesimal variation). 

The dynamics of the system is governed by the Lagrangian 

L = Lh + Lem + Lp + Lint (9) 

with 



LuigaP-, dxQap, dx^iQap) = \ Rds 

LEM{gai3,A^,daA^) = --^/\g\Ff"^Fp^ds 

LF{J^,9ap) = -V\9\ [P (1 + e(p))] ds 
Lint = qJ^Afj^ds 

where R is the scalar curvature of the metric, Fpo- are the curvature coefficients (field 
strength) of the connection A (gauge potential) and e{p) is a generic function of the scalar 

Iff T/i JV 

p = . /:^iii_ — ^ physically interpreted as the internal energy of the fluid. It gives rise to 

the pressure P = p^^. 

The gauge-natural bundle of parameters we introduce is then 

^(P) = (n^^p^yVLoriM) x^(p) {n^^p^TVCiV) x^(p) (7r^(p))*rM 

and the parametrization is such that for any section e = {6g^'^ ,6Ap,X^^) € T{E) and for 
any configuration a G r(C) we have 

<P5|/e>= {6gf,u{x'',gai3),SAp{x,Ac,),£xJ^) 
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where £xJ^ denotes the formal Lie derivative of J, i.e. £xJ^ = V,yJf^X'^ — J^V^X^^ + 
J^V yX" . In this way we are freely varying the gravitational and electromagnetic degrees 
of freedom, while variations of the fluid variables are taken to be tangent to the constraint 
S; in fact for every vector field X and for every admissible configuration J (i.e. such that 
^^J^' = holds) we have d^£xJ^ = 0. 

Let us remark that this parametrization is trivially gauge-natural; in fact it partially 
coincides with the identity and partially with the Lie derivative. 

Remark 33 We want to spend some words to stress that our parametrization is adapted to 
the constraint S, but in general not faithful to it; nevertheless this is exactly what Physics 
requires. Varying J according to our prescription is equivalent to drag J along a 1-parameter 
family of diffeomorphisms generated by X and this amounts to drag the integral curves of n'^ 
and to adjust the density p to keep J conserved (see [13]). If we consider a compact subset 
D C M and find a tangent vector X whose flow moves the integral curves of u'' also on 
the boundary of D {X\q£) ^ 0) but leaves their tangent vectors fixed {£xJ'^\dD = 0, as it 
is possible in some cases), we have that it produces a family of deformed J that are fixed 
on the boundary and that fulfill conservation of matter, but do not come from a vector field 
that vanishes on the boundary as it would be requested by faithfulness. Nevertheless, this is 
exactly what we want to do from the physical viewpoint: when we think of the congruence of 
curves that represent our fluid and we imagine to vary them without moving the boundary we 
do not want to move particles, not only tangent vectors. To support our choice to vary fields 
along our parametrization, we stress that also in the alternative approach of "abstract fluid 
particles" variations leave unchanged the particle identification on the boundary. Solutions of 
the Euler-Lagrange equations that we are going to find will thence be P5-critical sections of the 
parametrized constrained variational problem, without being necessarily vak-critical solutions 
of the variational problem with the constraint S given by dJ = 0. 

The first variation formulae arising from the freely varied Lagrangians Lh and Lem are 
the well known 



{6Lh \f < fs I j^e > > = ^{R^.u - \Rg^.u)5g^'' ds + V,(F-^) | j^ < P5 | jh > > (10) 

with 

(Ff^) I ji < Fs I j'e > > = ^ia'^^dpa - <5f/^))VA<55^'^ ds 
where the covariant derivative is relative to the Levi-Civita connection and 

{SLem \j^ <Fs\jh>) = - ^\{V^F^'nSA, ds - v^V^(F'^'^M,) ds 

(11) 



^4f^^)55^^d. 



with 

VW\Hif/^) ds = -2^^ = VW\{F,,F^.^ - -F^,F-Pg,y) ds. 



^ ^^ j^ = ^(p F/-i. 
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From the parametrized variation of Lp, by taking into account the definitions 



2^6 



fi = p{l+e{p)) and P = p — 



together with their consequences 



Pfp=f^ + P and pV.(^^)=V.P 



we have 



{6LF\<rs\j'e> 



dp 
dp 



u^ {VuJ^x" - rvuXf + j^VuX") ds. 



Integrating by parts, we get the first variation formula 

{6Lf \<^s\J^e>) = -^\ [{u-^u^ + 5^^)V,P + {p + P) w^V.uJ X^' ds+ 

H + P 



P 

Is] 



u^X^iV^J^') + vl^V. [{p + PW^ + u^u'')X^'] ds+ 
HpSg'^'^ds 



(12) 



with 



hI^^ ds = 2^-^ = vl^l [POfiu - (m + P)K'^u] ds. 



Let us compute the first variation formula for the interaction term L^nt)'- 



[5L 



int 



< 



\j £> 






dJ^" 



dA, 



qA^{V^J^X'' - J^V^X^" + Ji'V^X") ds + qrSA^ds 



(13) 



= qrSA^ + qJTf.^X'' ds + gV^(yl^ J^X^ - A^rX^") ds 
The complete set of the equations of motion for the Lagrangian (9) is thence the following 
(S/^J^' = (constraint) 



R 



ii?<,,. = .«*^) + <)) 



~g\V,F>^'' = qr 



[ {u^u'' + jpV.P + {li + P) u'^V.u^ - Jf^Ff,, = 0. 

We are now going to prove that the Lagrangian (9) is variationally gauge-natural when 
restricted to the prolongation of the constraint submanifold 5 C J^^{P) identified by 
equation V^J'^ = 0. 

If (x^, 9) are coordinates on the structure bundle P, an infinitesimal generator of principal 
automorphisms of P is a vector field Y = y^(x°)9^ + y{x°')-^ G T{TP) projecting onto 
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X = TttY = Y^d^. The Lie derivative of a section a of the gauge-natural configuration 
bundle Lor{M) x C{V) x K^T*M with respect to Y is 

n r\ O 

£yCJ = {£x9a(3)T^— + {£yA,)^:^ + {£xJ^) 



Ja/3 



dA, 



5J^ 



with coefficients 



£xgai3 = v^y^ + V^F„ 

£yA^ = YPFp„ + V^{y + ApYP) 

£yF^P = YPVpF^p + VoX'Fpp + VpYPF^p 

x^ji^ = Vyj^Y" - rv^Y^ + Jpy^Y". 

Thus the configuration bundle is gauge-natural with order (1, 1). 

A gauge-natural morphism J : J^d — > (J^IGA(P))* ^jicSl-P) such that we have 



<P5 I j' < i\j^Y > > o /(J = £Ya 
can be defined by < J|i^l^ > o j^a = {£xgai3, £YAa,X). 

Let us now consider the Lie derivative of the morphism L with respect to Y. We find: 

dL 



(14) 



TT2£yL =- £xgaP + 777^ d\£xgali + 

Odxgaf3 



-d\uj£xgai3+ 



dgaf3 ddxgap ddx^gafs 

r)T BT BT 

Analyzing this term by term we have that the identities £yLh = 0, £yLem = are well 
known to be satisfied (see [7]) and they can be recasted using the first variation formulae 
(10) and (11) together with the identity (14) into 



with 



and into 



Off as a 



2k 



2k 



-(R, 



1 



111/ 



Rgnu)£Yg'"' = V,f_ 



-R"x 

2 ^ 



i?<5? Y^ + 



H 



«_^«(7<5j))v^^y^ 



dSa 



(15) 



\g\{V^F^^n£YA. - ^H\f/'^£xg^'' = V^£% 



■-EM 



with 



£em d^f 



gua jjiEM) yn ^ i.<^/.v,,(y + ApYP) 



ds. 



(16) 



For the Lagrangian Lp a straightforward calculation enables us to verify that the identity 
£yLf = holds too, and that it can be split (thanks to the first variation formula (12) 
together with identity (14)) into the following: 



\g\ [{upu'^ + ^"iD^uP + (m + P) u^v.u^] XP+ 



1-1 + p 



u^X^iVpJP) 



^-HJ£^£xg^'' = Va£^ 



(17) 
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with 

^^ dsa = - vl5l (5"^ H^ Y^) ds^. 

What remains is just to show that the Lagrangian Lint is variationally gauge-natural when 
restricted to the constraint submanifold. In fact we have 



7r2 £yL' 



int 



^^'^'£xJ^ + ^£yA. - V^{Y-Unt) 



ds 



dJt" dAy 

= q [A^(V,J^y" - rV,Y^' + J'^V.y^) + r[YPFp, + V,(y + ApYP)] + 

- A^ry^Y"" - Y'^ry^A^ - Y'^A^Vo^r] ds = 

= q [rV,y] ds = q [V,(y J") - yV^J"] ds 

where the first addendum is a (local) divergence, while the second one vanishes on the 
constraint. Thanks to the first variation (13) together with identity (14), the previous 
formula can be recasted into the identity 

qr£YA, + qJ^Fp.X'^ = VaS^nt (18) 

with 

£rnt = -q{y + ApYnr 

Summing up identities (15), (16), (17) and (18) we find that the sum of the left hand sides 
vanishes if composed with a P^-critical section, while the current 

£ dSa = {£h + £eM + £f ~^ ^int) ^^a 

is conserved. 



If in particular we compute the current £" dsa relative to the horizontal part y(hor) 
Y^dfj_ — ApYP-^ of an infinitesimal generator of automorphisms of P, what we get is 

^fhor) ds^ = £% dSa - V¥\ [5"^ (HP + HP)Y''] dSa. 

and this is the conserved current related with the energy-momentum of the system. 

'{v) = {y + Ap} 



Moreover, the Nother current relative to the vertical part Yiy\ = {y + A^y^)^ is 



£?V) = - VM ^"''V^(2/ + ApYP) -{y + A^Yf') q,r. 



Integrating by part the covariant derivative of {y + ApY^) we can split this vertical current 
according to 



^iV) = (\^V^F"^ - q,n{y + A^Y^ - V4^^F"^(y + ApY^)] 

where the first summand vanishes on-shell while the argument of the divergence is called 
the superpotential U = —^y'\g\F°'P'{y+ApYP)dsa^ and it is a 2-form that is closed even off- 
shell. We stress that the interaction Lagrangian does not contribute to the electromagnetic 
superpotential. 
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Conclusions 

Motivated by many physical and mathematical examples, we have studied variational 
problems with parametrized variations. In particular, starting from the results of [8] we 
have worked on a twofold generalization: by one side we studied parametrized variational 
problems, no matter from where the parametrization comes from, and we have recovered 
and ever generalized known results. In a second moment we also focused on the applica- 
tion to vakonomic constraints, for which we have defined and compared the two different 
concepts of vak-criticality and P-criticality. A detailed study of nonholonomic field theory 
could also be based on the same framework and it shall certainly be performed in the 
future developments. Hopefully our unified language will also help to understand how 
constraints on the derivatives of fields have to be handled in the case of field theories. 
For what concerns the second direction of generalization, we have formulated a variational 
theory for gauge-natural parametrized field theories including conserved currents that 
provides, in the case of relativistic hydrodynamics of a charged fluid, a non-conventional 
procedure to define field equations, Nother currents and superpotentials and that turns 
out to be computationally much easier then the standard one. A last observation that 
deserves to be further investigated is that a P-Euler Lagrange equations arising from 
the parametrized variational problems is usually not variational in the standard sense. 
Nevertheless we have a very efficient machinery to study symmetries and conservation 
laws. Do this apply to any known non-variational equations? 
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A Variational morphisms and their splittings 

In this Appendix we report the theory of variational morphisms. The theory is modeled 
on [7] and [10], where complete proofs are given. To adapt them to the constrained case 
we have to modify a bit the definitions. In particular, the modified definitions of formal 
connection and variational morphism seem to appear here for the first time. 

A.l Variational morphisms 

Variational morphisms are an abstract model of the integrands appearing in global vari- 
ational calculus. They were introduced (see [6] for a discussion and some bibliography 
on alternative approaches) to provide a general framework in which we can implement 
two specific algorithmic procedures based on integration by parts. The first one explic- 
itly construct a splitting of such integrands into a volume and a boundary covariant part 
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that generalizes the first variation formula; the second provides superpotentials from the 
Nother currents. 

Definition 34 Let C — > M be a configuration bundle over a m dimensional base and 
E — ^ C a vector bundle over C. Let k > h and n be natural numbers (zero included). 
A variational E-morphism is a vertical (projecting onto the identity) bundle morphism 

M : J'^C — > {f'E)*(^jhcT^*K"'"'''T*{M). 

The minimal k is called the order of M, h is the rank, while m — n \s the degree and n the 
codegree. 

A. 2 Operations with (variational) morphisms 

Definition 35 Given a vertical fibered morphism F = (id,<l>) between two bundles B and D 
over the same base M we call s-order jet prolongation (or s-order formal derivative) of <I> the 
unique fibered morphism j^'F = (id,j''$) that makes commutative the following diagram 



f<y 



^J'D 




fi^oa) 



When one prolongs a morphism projecting onto the identity from a jet bundle to another 
bundle, e.g. F : J^B — > D, one gets a morphism j^¥ : J^{JPB) — > J^D. The bundle 
jp+^Q is naturally embedded in J^{J^B) by a canonical inclusion i. We are often more 
interested into the map j^¥ o i rather then in the prolongation itself. We call j^¥ o i the 
first order holonomic prolongation of F. In the sequel, when we mention the prolongation 
of a morphism starting from a jet bundle we will always mean the holonomic one. 

Remark 36 We warn the reader that by j'' we denote three different concepts of prolonga- 
tion: prolongation of a section, of a vertical morphism and of a vector field. They are strongly 
linked, though not identical concepts. For a discussion of this link we refer the reader to [7], 
but we just remark that every vector field X on a bundle B can be seen as a vertical morphism 
X : B — > TB, and if we compute the first jet prolongation of this morphism we get another 
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morphism j'^X : J^B — > J^TB that is not a vector field on J^B as it should be if thought 

as the prolongation of a vector field. 

If (x^, y]j) are coordinates on J^B {v is a multiindex of length < jz^l < h), and (x'^, y^, Up , f?) 

coordinates on J^TB, the coordinate expression of the prolonged morphism j^X \sj^X{xf^, yl) 

{xl',yl,dpv^',dpv''). 

Luckily enough, for all integers h, a global canonical map r^ : J^TB — > TJ^B can be 

defined (see [19, 12]), whose coordinate expression is rh{x'^,yl,Vp,Vp) = {x'^,yl,v'^,dpv''' — 

yZnV^-^i dp-aV^) where a is any possible multiindex of length < lal < |p| = h and 1„ is 

a multiindex with 1 entry in position /i and zero elsewhere. The ordinary prolongation j^''X of 

the vector field X as a vector field on j'''B can be recovered just by composing rh with the 

prolonged morphism j^ . 

Definition 37 Let B — > M be a bundle and E a vector bundle over the same base, let M 
be a fibered morphism M : B — > (E)* and P another bundle morphism P : B — > E. 
We define formal contraction of M on P the unique variational morphism 

<M|P>:S — > M xR 

such that 

Vp G r{B), <M|P>op=<Mop|Pop> 

where < | > in the right-hand side is the usual fiberwise contraction between a section of a 
vector bundle and a section of its dual. 

Definition 38 Let M : j'^C — > A'"-"r*(M) be a variational morphism. We define the 
divergence of M to be the unique variational morphism Div(M) : J^+^C — > A™~"+^T*(M) 
such that for any section a € r(C) the following holds 

Div(M) oj^+V := d(MojV). 

A. 3 Formal connections and fibered connections 

The following definitions generalize the ones given in [10] to the settings of constrained 
variational calculus. Let us consider a composite projection E — *• B — > M and the 
fibered morphism (section) 



> E 




O O" 



with E vector bundle on B. 

Let {x^, y*) be a coordinate system on B and (x^, y*, u ) on E (with respect to a fiberwise 

local basis {eyi}) a local representation of F is u^(x^,y*)eyi, while a local representation 
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of its first jet prolongation (as a fibered morphism) j^F is v {x^,y^)eA + d^v e^, where 

{e^ is a base of the fiber of J^E — > E and df^v = d^v + diV y'' is called the formal 

derivative of F. 

The transformation rules of the formal derivatives under a fibered change of coordinates 

are not tensorial. 

Consider a set {T^^{x'',y\yl,)}A,Be{i,- ,rankE}, ,ie{i,- ,dimM} of local coefficients that fit 

into the expression 

giving rise to an object at the left hand side of the equalities that under a change of fibered 
coordinates 



yy = dy + Viy (19) 






xt" = x" (x") 

, v^' = Mf (x^^, y')v^ Mf G GL(rank^, M) 

transforms according to the rule 

V^.t;^' = {J-X'^^yui' (20) 

that are nothing but the transition functions of the vector bundle E (Sib T*M — > B. 
The resulting transformation law for the coefficients is 

r^V = j;,(r^^Mi'M|, - d.Mi'ME,) (21) 

It's easy to verify that relations (21) form a cocycle over the manifold J^B, thus they 
define the transition functions of an affine bundle FC{E) — > J^B modeled on the vector 
bundle E* ®b E ®b T*M. 

Definition 39 A formal connection on the composite fiber bundle E — > B — > M is a 
global section of the affine bundle FC{E). 

Let us consider the bundle L(£') — > B of linear frames of the vector bundle E — > B. 
Let r = rankE'. L(£^) is a principal GL{r,M.) bundle and the free natural right action 
r : L(^) X GL(r,]R) — > l.{E) can be defined on it. 

Let us call "Lm^E) the bundle L(£') — *• M obtained by the obvious composition of projec- 
tions. It is no longer principal, nevertheless the action a induces on it a free right action 
ttM-, that can be lifted to a right action j^om on its fist jet extension J^L,m{E). The 



action j^om is still free and admits a quotient manifold qjjI.^) that has an affine bundle 






structure over the base J^B. The two affine bundles FC{E) and qijIS^\ are isomorphic 
and this can be proved directly by showing that they share the same transition functions, 
providing a more intrinsic characterization of the bundle of formal connections. 

Consider now the bundle J^Eji^ — > J^B; the presence of a formal connection on E, 
allows us to define the global vector bundle morphism V by the following commutative 
diagram 
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J^Ejib^E(E)bT*M 



1 
J^B ^^^ B 



whose local coordinate expression is (19). 

Definition 40 Let a be a section of J^Eji^- We call Vocr the formal covariant derivative 
of a. 



■^E 



Definition 41 A fibered connection on E — *■ S is a couple (7,r), where 7 is a linear 
connection on M and F is a formal connection on E. 



Given a fibered connection on E we can extend the morphism V to J^{Ejig ® Tq{M)) 
where Tq(M) is the algebra of p-times covariant and g-times contravariant tensors on M 
with the obvious tensorization procedure, getting a fibered morphism 

JHEjib ® T^{M))^^E (8)B T^+i{M) 



J^B — 

that we denote by the same symbol V. 



^Bo 



-> B 



Definition 42 Let ct be a section of J^ {E ji B®Tq {M)) , let us call Voa the formal covariant 
derivative of a. 



Definition 43 The formal covariant derivative of a morphism F : B — > E is the unique 
fibered morphism VF : J^B — > E®b T*M such that 

V/9 G T{B), VF o j^p = V o ji(F o p). 

A formal connection on a composite fiber bundle E — > B — > M can be constructed from 
a linear connection A on the vector bundle E — > B and the first jet prolongation j^a of 
a section of the bundle B — > M. Let (x'^,y",u ) be a fibered coordinate system on E; 
the connection A has the coordinate representation 



A = dx^ ® {^ - A?„.v^ ^ 



e^,--D,^ 9^) + ^y^^&-^g.-^5^ 



d_ 



)• 



The coefficients r£ = A^ + A^-y^^ transform according to (21) and thus they glue 
together into a global formal connection on E — > B — > M. 



29 



A. 4 Local expression and reduction of a variational morphism with re- 
spect to a fibered connection 

Fixed a fiberwise basis s eyi(x), e^^(x), • • • ,e^^" ''(x)> of 

L J Ag{1,--- ,p}, iGJl,--- ,/i}, Aig{l,--,m} 

J!I:E and its dual \e^{x),ej^ (x),--- ,e^...;^ (x) > in {J^E)*, the local expression of a 
variational morphism M : j'^C — > {J^E)* ® A"'-'^T*{M) is the following 



n\ 



where the indices //i • • • ^n are skew-symmetric while the Xi- ■ ■ Xh (if any) are symmetric 
and where the coefficients are functions of the point j^p G J C. 

Given a fibered connection {iguJ^Bu^ ^^ ^^^ change basis on J^E in a way that a holo- 
nomic section j X € T{J E) in the new basis (called basis of symmetrized covariant 
derivatives) has components 

vf^ = v,y{x) = d,y + ri,y 



<...A„=V(A,---V,„)t;^(x). 

Accordingly, using the dual basis on {J E)* we can define a new coordinate system 
{x°',v^2'^'',v^/'^''^\v^^'"^'"^^'^'') on {f'Ef (g)jhc^"'~"'T*{M) such that the local ex- 
pression of M is 



^y^A ^ -^""A "=^Ai -r • ■ ■ -r ^A ^Ai-A^y '^"•^m-Mn- 

We remark that in the previous expression, as in formula (22), the coefficients have to be 
meant as functions of the point j^p G J C, and the indices pi- ■ ■ pn are skew-symmetric 
while the Ai • • • A/i (if any) are symmetric. 
The explicit change of basis and coordinates in the case /i = 1 is 

1 pA _ A.yA B ^'^^ \ A _ ^A 
y^\ — ^X + ^ BX^ [^A — ^A 

The introduction of a fibered connection can select a "preferred" kind of variational mor- 
phism having nice properties when written in the symmetrized covariant derivatives coor- 
dinate system. 

Definition 44 Let M : j''C — > (J^E)* (^jhc A"'~''T*{M) be a variational morphism and 
M/ its term of rank < I < h. The term M/ is said to be reduced with respect to the fibered 
connection (7fl„,r4 ) if, written in the symmetrized covariant derivatives coordinate system, 

it has the coefficients symmetric in the first n + 1 upper indices, i. e. if {j^^^'"'^'^ ^'" ' = 0. 
The variational morphism M is reduced if each one of its terms is reduced. 
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Notice that the term of rank zero is always reduced with respect to any connection, 
having always antisymmetric upper indices and that in the case of Mechanics (m = 1) all 
variational morphisms are reduced with respect to any connection. 

Now we are ready to state the most important results of the theory, that justify the name 
"variational" for the morphism we have introduced, providing a sort of "first variation 
formula" from any variational morphism. The proofs of the next two Theorems can be 
found in [7], but as a consequence of the minor modifications we introduced here one has 
to reinterpret the symbols according to the previous definitions. 

Theorem 45 -"Splitting lemma"- Let M : J^C — > {J^Ef ®jhc N^T*M be a global 
variational morphism with codegree n = and let /i, /c G N (zero included). Chosen a fibered 
connection (7, F) on E there exist a unique pair V, B of variational morphisms 

V = V(M) : J^^^C — > E* ®c A"'T*M 

n = ]B(M,7) : J^^^~^C — > {J^^^Ef ®jh-ic A"'^^T*M 

reduced with respect to (7, T) such that Vy € r(£') the following holds true 

< M I j^V >=< V I F > + Div < B I /"V > . (23) 

The variational morphism V is called the volume part of M and B its boundary part. 

Theorem 46 -"Reduction lemma"- Let M : J^C — > {f'E)* ®jhc A'^~"-T*M be a 
global variational morphism of codegree n > 1 and let /i, A: € N (zero included). Chosen a 
fibered connection (7, F) on E there exist a unique pair V, B of variational morphisms 

V = V(M) : j'^^^C — > (J^E)* ®jhc K"'~''T*M 

B = B(M,7) : Z'+^'-^C — > {f'~^E)* ®jh-i A"*-"-ir*M 

reduced with respect to (7, T) such that \/V € T{E) the following holds true 

< M I /y >=< V I j^V > + Div < B j i''" V > . (24) 

The variational morphism V is called the volume part of M and B its boundary part. 

Constructive proofs of these Theorems are given in [7] where we can find a collection of 
older results; they are carried on by induction on the rank of the morphism. If the rank 
is 1 then one just needs to integrate by parts and no fibered connection is needed. At any 
rank they provide an algorithm out of which we can perform explicitly the splitting or the 
reduction. 

In both the Theorems the volume part is uniquely defined even if we do not require it to be 
reduced, and finally it does not depend on the fibered connection. The boundary parts, on 
the contrary, are defined modulo a divergenceless term, but the condition of being reduced 
with respect to any specific fixed fibered connection determines them uniquely (if the rank 
is 1 the boundary part is stiill unique and reduced with respect to any fibered connection). 
Moreover it was proved [15] that the boundary parts depend in fact only on the connection 
on the base, while the formal connection on E reveals itself to be just an intermediate 
object useful to explicitly construct the splitting, but in the end unessential. 
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B Gauge natural bundles 



Definition 47 Let G be a Lie group. *P"(G) will denote the category of principal G-bundles 
whose base manifold is n-dimensional and whose maps are G-bundle maps projecting onto 
local diffeomorphism between the two bases. 

Definition 48 A covariant functor 53 from the category *P"(G) to the category of bundles 
and bundle morphisms is a gauge-natural functor '\f 

a) for all principal G-bundle P — > M, 5S(P) is a bundle with structure group G having 
the same base M as P; 

b) any principal morphism <I> : P — > P' projecting onto (/) : M — > M induces a fibered 
morphism *B($) : *B(P) — > *B(P') also projecting over 0; 

c) for any open U C M the inclusion morphism i : p^^{U) — > P is mapped into the 
inclusion morphism <B(i) : ^{p~'^{U)) — > 53(P). 

Definition 49 A gauge-natural bundle GN — > M having the principal bundle P — > M 
as structure bundle is the image of P trough a gauge-natural functor iB. 

Examples of gauge-natural bundles are all natural bundles (trivial case), the bundle of 
principal connections C{P) on a principle bundle P, the bundle IGA(P) of its Infinitesimal 
Generators of principal Automorphisms, the bundle of spin frames and others (see [7]). 
Notice that jet prolongations of a gauge-natural bundle are gauge-natural with the same 
structure bundle, and that every composite fibration Z — > Y — > X where y is a gauge- 
natural bundle with structure bundle P and iiz is natural makes vry o vr^ a gauge-natural 
bundle with structure bundle P. 

We want to stress that for our purposes, the most relevant property of a gauge-natural 
bundle 53 (-P) is the possibility of lifting local principal automorphisms of the structure 
G-bundle P to local fibered morphisms of 53 (P). We will call ^^ gauge transformation^^ on 
03 (P) the gauge-natural lift of a local principal automorphism of P and ^^ gauge groups the 
structure group G. 

The same, via the tangent map, can be done for infinitesimal generators: an infinitesimal 
generator of principal automorphisms can be lifted to a vector field on 53 (P) that will be 
called its ^''gauge-natural liff . 

Let us consider a 1-parameter family of infinitesimal generators of principal automor- 
phisms {^ s = ('^sjV's)}- It can be proved (see [7]) that any infinitesimal generator 
H G r(rG) = ^^sl^^o is a right invariant projectable vector field on P and that we can 
explicitly construct the bundle IGA(P) — > M whose sections are in one to one correspon- 
dence with the infinitesimal generators of such families and with right invariant sections 
of TP. Infinitesimal generators of vertical principal automorphisms {\I'<j = (id, f/'s)} can 
be identified with sections of a subbundle ^jr- ^-> IGA(P) that is associated to P — > M 
by means of the adjoint action of the Lie algebra g of the gauge group G. 
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By means of a trivialization of P we can locally define a set {pA}AG{i--dimG| of right 
invariant vector fields on P that at every point form a basis of VP. Let (x^, y") be a fibered 
system of coordinates on P; every infinitesimal generator of principal automorphisms H 
can be locally written as a right invariant vector field 

on P. 



Definition 50 We say that a gauge-natural bundle 5S(P) is of (r, s)-order [s < r) if the 
gauge-natural lift is a morphism 

L : «B(P) — > (jnGA(P))* 02J r«B(P) 

such that if H is vertical (i.e. it is a section of ^ ^^ IGA(P)), the vector field E<q =< 
L I j^'H > on 5S(P) is vertical, too, and depends in fact on the derivatives of E only up to the 
order s. 



Definition 51 Let ^ = {4>,iIj) be a local principal automorphism of P and *B(^) = ((^j^rg) 
be its gauge-natural lift to the gauge-natural bundle 5S(P). Let cr be a section of *B(P); we 
define its push-forward along ^ to be the section ^!*a = -0*3 o a o (f)~^ . 

Definition 52 Let H be the infinitesimal generator of the 1-parameter family of local principal 
automorphisms {^ s = {(f)s,iljs)}; we define the Lie derivative of a section a of 5S(P) along 
H to be the unique section of the pull back bundle £=0" € r((T*F*B(P)) such that 

Remark 53 We can easily prove that if Hsg is the gauge-natural lift of H on 5S(P), and H is 
the projection of H on M , the Lie derivative of a section can be computed trough the following 
formula 

£■^0 = Ta o ttH — Hgs. 

In coordinates, if {x^,y^) are fibered coordinates on 53(P) while the lift can be written as 

It also turns out that the Lie derivative operator on sections of the gauge-natural bundle 53(P) 
of order (r, s) (where r > s) can be interpreted as a fibered morphism 

B£ : J^^{P) — > {r\GA{P))* V^{P) 

such that for all sections a G r(!B(P)) and for all infinitesimal generator of principal automor- 
phism E e r(IGA(P)) one has < B^ | J^'H > o jV = £e(t. 
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Notice that if H is vertical (i.e. it is a section of ^ ^^ IGA(P)) then £=c7 depends in 
fact on the derivatives of H only up to the order s. 

Remark 54 When in differential geometry one computes the Lie derivative of a vector field 
on a manifold M one gets another vector field and so happens for differential forms and tensors. 
According to our definition, on the contrary, the Lie derivative of a vector field is a section of 
the vertical bundle VTM. To reconcile this apparent contradiction we have to remark that 
whenever *B is a vector bundle with base M then its vertical V^ is isomorphic (see [16]) to 
the fibered product OS xj^f *B and one has a global vertical fibered morphism 7r2 : V^ — > OS 
that realizes the projection onto the second fibered factor. If {x'^,y^,v"') is a local coordinate 
system on y*B then the map 7T2 acts as follows: 7r2 : {x^^,y"-,v°') 1 — > (x^,u"). If cj is a 
section of a gauge-natural vector bundle then 'K2£'e(t would be what traditionally is meant 
as its Lie derivative, the coordinate computation being identical, but with a slightly different 
interpretation. 

Let P — > M be a principal bundle; let moreover U d M he, open, [x^,g'^) be a coordi- 
nate system on p~^{U) C P and {pa{^^ ^ 9'^)} A=i-- dim G ^ fiberwise right invariant basis 
of vertical vector fields. An infinitesimal generator of principal automorphisms has the 
local form E{x'^,g'^) = ^{x)'^dfj_ + £,^{x)pA{x^,g°') (see [7]). Given a principal connection 
u{x'^, g"-) = dx^®{df^—u>^{x)pA{x^,g"')) we can also split H into a vertical and a horizontal 
part according to 

S = eidf. - uj^Pa) © (<e + ^^)PA = H(H) e H(y). 
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